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We calculate the tunneling current for a bilayer quantum Hall system in the interlayer incoherent 
regime. In order to capture the strong correlation effects we model the layers as two Wigner crystals 
coupled through interlayer Coulomb interactions, treated in the continuum limit. By generalizing 
previous work by Johansson and Kinaret (JK), we are able to study the effect of the low energy 
out-of-phase magnetophonon modes on the electron "shake-up" which occurs during a tunneling 
event. We find the tunneling current peak value to scale with the magnetic field as found by JK; 
however, we find a different scaling of the peak value with the interlayer separation, which agrees 
with the measurements by J. P. Eisenstein et al, Phys. Rev. Lett. 74, 1419 (1995). 
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I. INTRODUCTION 



As a consequence of interlayer interactions, bilayer two 
dimensional electron systems can exhibit novel behaviors 
in the high magnetic field (quantum Hall) regime. The 
virtue of the bilayer system is that the relative strength of 
the intralayer to the interlayer interactions can be varied; 
their relative strengths are characterized by the dimen- 
sionless quantity d/l, with d the interlayer separation and 
I the magnetic length. For example, a balanced bilayer 
system (i.e, with equal densities) at total filling factor 
vt = 1/2 (each layer in a v — 1/4 state) or vt = 1 
(each layer in a v = 1/2 state) exhibits quantum Hall 
plateaus which would have been forbidden for the single 
layer systemic A spectacular manifestation of interlayer 
interactions for vt = 1 is the onset of spontaneous inter- 
layer coherence, which occurs for closely spaced layers 
(with d/l < 1.8). The interlayer coherent regime is char- 
acterized by a prominent zero bias peak in the tunneling 
conductance, 3 which disperses with the application of an 
in-plane magnetic field. 4 Recent transport studies pro- 
vide compelling evidence that the coherent state can be 
understood as a bilayer excitonic condensates^ 

Here we will focus instead on the tunneling properties 
in the incoherent regime. Even in the absence of inter- 
layer coherence there are still important effects of the 
interlayer interactions, in particular the Coulomb barrier 
peak 3 ^ which is observed in the tunneling current-voltage 
characteristic (I-V). This behavior is a consequence of 
both the interactions and the electron dynamics in high 
magnetic fields. A semiclassical picture is that the elec- 
trons in the accepting layer accommodate a tunneling 
electron by forming a correlation hole. However, the cy- 
clotron motion of the electrons in a high magnetic field 
frustrates the formation of this hole (rather than mov- 
ing radially outward to form the hole the electrons in 
the accepting layer precess, due to the Lorentz force), 
and the tunneling can only proceed if the tunneling elec- 
trons are externally provided with sufficient energy to 
overcome this "Coulomb barrier". The result is a peak 
in the tunneling current at a bias equal to the typical 
Coulomb energy, such that eV w e/Aireao, where ao is 



the average spacing between electrons. The width and 
shape of the I-V characteristic will depend upon the in- 
terplay of intralayer versus interlayer interactions (i.e. 
on the parameter d/l), and is the subject of this paper. 
The semiclassical picture described above is essentially a 
dynamical version of a static exciton picture of the tun- 
neling process, first discussed by Eisenstein et al.^, in 
which the tunneling electron and the hole that it leaves 
behind combine to form an exciton; hints of the exci- 
ton's existence emerge in shifts of the I-V characteris- 
tics with interlayer spacing (the Coulomb barrier peak is 
"red-shifted" to a lower potential by virtue of the inter- 
layer Coulomb energy —e 2 /ed). 

To study the tunneling current we adopt a model pro- 
posed by Johansson and Kinaret- (JK), in which the two 
dimensional electron system in each layer is assumed to 
have formed a Wigner crystal (WC). Of course, for the 
experimentally realized filling factors the layers are in a 
liquid phase; the JK model is only intended to capture 
the important short distance correlations among the elec- 
trons, and the magnetophonon modes of the WC provide 
a reservoir of low energy excitations that can dissipate the 
energy of the tunneling electron. In the approach of JK 
the two layers were assumed to be noninteracting; we de- 
part from their approach by incorporating the neglected 
interlayer interactions. These interactions qualitatively 
change the collective mode structure of the bilayer system 
by introducing a gap in the out of phase magnetophonon 
mode and by affecting dramatically the I-V character- 
istics of the bilayer system since, as we mentioned ear- 
lier, they "red-shift" the peak bias values by an amount 
proportional to the excitonic energy e 2 /ed. The tunnel- 
ing electron couples to the out-of-phase magnetophonon 
mode, and the electronic "shake-up" which accompanies 
the tunneling event depends upon the behavior of this 
mode. 

This paper is organized as follows. In Sec. II we re- 
view the continuum theory of the single layer WC. In 
Sec. Ill we generalize this to the bilayer case, including 
both long range and short range correlations between the 
two layers, and find the collective modes for the bilayer 
WC. We also determine the coupling between a tunnel- 
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ing electron and the collective modes of the bilayer WC, 
which can be written as a type of spin-boson model^S 
Section IV is devoted to the calculation of the tunnel- 
ing current, which is expressed in terms of a correlation 
function that is determined by solving a linear integral 
equation both numerically and with an approximate an- 
alytical method. The solution allows us to qualitatively 
reproduce the tunneling current curves and obtain a good 
theoretical estimate for the peak bias values. It also re- 
produces the expected peak bias voltage dependence of 
the I-V curve with interlayer distance. Many of the de- 
tails of our calculations are provided in the appendices. 

II. SINGLE LAYER THEORY 

We will start by introducing the single layer model, 
in which the 2-D electron system is assumed to be in a 
Wigner crystal state, modeled as an elastic medium with 



a momentum cut-off of go = 2y / 7mo, where no is the sin- 
gle layer density (the cut-off is determined by imposing 
conservation of the total number of states). Of course, for 
the filling factors relevant to the tunneling experiments 
this is not the ground state of the true system. This sim- 
plified approach to the highly correlated liquid-like state 
of the electron system takes care of the short range corre- 
lations among electrons and creates an abundance of low 
energy collective excitations, which are needed to relax 
the local defects created by tunneling events. In order 
to describe the system as an elastic medium we need to 
introduce the two Lame coefficients (the compressional 
and shear moduli) A, fi which are associated with the en- 
ergy cost of deformations in the systemiii The effect of 
the perpendicular magnetic field is introduced through 
the vector potential A. The Lagrangian that describes 
the dynamics of the 2-D electron system in the presence 
of the perpendicular magnetic field and the intralayer 
Coulomb interaction is given by 
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where e is the dielectric constant of the host materials 
(GaAs in this case) and u is the electron lattice displace- 
ment field. Notice that for the intralayer Coulomb inter- 
action term we have used the linear continuous approx- 
imation for charge density fluctuations Sn/riQ — — V • u 
(this is correct in the absence of vacancies and inter- 
stitials). We choose to work in the symmetric gauge 
A(u) = (—Buy/2 , Bu x /2 , 0) where B is the applied 
magnetic field. In the absence of B the eigenmodes of 
the elastic medium can be labeled as either longitudinal 
or transverse; the magnetic field couples those two modes 
together and one has to diagonalize the Hamiltonian of 
the system anew. Nevertheless it is useful to decompose 
the displacement field u into transverse and longitudinal 
components ut and ul- After Fourier-transforming, the 
above Lagrangian becomes 



and oj c — eB /m is the cyclotron frequency. If we compare 
these modes with the results of Bonsall and Maradudin 15 
for a WC with pure Coulomb interactions, we see that 
cl is a higher order correction to the magnetoplasmon 
contribution (the second term on the right hand side of 
Eq. |3J), and can in most cases be neglected. 

We show in Appendix A that the new eigenfrequencies 
for the collective modes of the single layer electron system 
are given by 
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as the longitudinal and transverse zero magnetic field 
eigenfrequencies with 



cl 



\/{\ + 2n)/mn , 



mno, 



In the zero magnetic field limit (lu c — 0) the above modes 
decouple into pure longitudinal and transverse modes, 
while for q — * 0, the magnetoplasmon mode uj + — * to c , 
in accordance with Kohn's theorem^ For lo c 3> lul, ljt, 
the eigenmode frequencies can be expanded as 



uJ7 



LJ C 



2uj, 
lol^t 



+ OK7 3 ). 



The magnetophonon mode w_ 



q 



3/2 



(G) 
(7) 



at long wave- 



(4) lengths, in agreement with Bonsall and Maradudin 
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III. BILAYER THEORY 



IV. TUNNELING CURRENT 



We now consider a bilayer system composed of two 
parallel single layers A and B, separated by a distance d, 
with the magnetic field normal to the layers. In a contin- 
uum model, the interlayer Coulomb interactions include 
a long range term which accounts for interactions be- 
tween density fluctuations in the two layers, and a short 
range term which accounts for the commensuration en- 
ergy between the WCs in the two layers. The latter term 
imposes an energy penalty for out-of-phase fluctuations 
of the electron densities in the two layers, proportional 
to (ua — ub) 2 - The Lagrangian for the bilayer system is 
then 



L = L A +Li 



n 
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IK 
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e 2 [V-u^(r)][V'-u B (r')] 



4ne^/(x - x') 2 + (y - y') 2 + d 2 



,(8) 



where La, Lb are the individual single layer Lagrangians 
given by Eq. J2J, and K is a "spring constant" for the 
out-of-phase fluctuations. When written in terms of the 
in-phase v = (u B + ua)/2 and out-of-phase u = u B — u A 
field displacement modes, the Lagrangian decouples into 
independent in-phase and out-of-phase terms, so that 
there are two in-phase and two out-of-phase eigenmodes; 
the details are included in Appendix B. The in-phase 
modes correspond to a uniform translation of both layers, 
and are unimportant for the tunneling processes which 
are the topic of this paper. The dispersion of the out-of- 
phase modes is similar to Eq. (|SJ) and given by 



2 



K 2 + 



n 2 T + n 2 L f - Att 2 L VL 2 T 



"I 1/2 



(9) 



where 



2 2 



2K 



2K 

mno ' 



e 2 n 
2me 



<?(1 



(10) 
(11) 



Both of these modes are gapped at long wavelengths due 
to the short range interlayer coupling K. In the high 
magnetic field limit we obtain for the magnetoplasmon 
and magnetophonon 
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So far we have developed an elastic theory for the elec- 
tronic charge fluctuations in the bilayer system, whose 
low energy excitations can provide an energy dissipating 
mechanism necessary for the system to relax any charge 
defect created by tunneling events. What remains to be 
done is to couple these collective modes to the tunneling 
electrons in the system. Following JK we assume that 
the tunneling events are sufficiently infrequent that we 
may treat them as occurring independently, and hence 
develop a model of a single tunneling electron which is 
coupled to the WCs in each layer. Such a model can in- 
corporate the effects of interlayer interactions (through 
the bilayer collective modes), but not interlayer coher- 
ence. The resulting coupling of a tunneling electron with 
the magnetophonons (the bilayer charge density fluctu- 
ations) is a standard electron— phonon coupling and is 
further discussed in Appendix C. The final form of the 
Hamiltonian is 



Hbath + Hrp + Hrp , 



(14) 



where 
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(15) 
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In the above M s a(b) are the electron-magnetophonon 
couplings and £a(b) the corresponding Madelung ener- 
gies for the two Wigner crystal lattices, c and <y are the 
tunneling electron annihilation and creation operators, 
a s and o| are the annihilation and creation operators for 
the collective modes, and T is the interlayer tunneling 
matrix element^ The collective modes that couple to 
a tunneling electron provide a mechanism to generate a 
"shake-up" that will relax the defect it has produced and 
dissipate the extra bias energy the tunneling electron has 
acquired. 

The tunneling current related to the above Hamilto- 
nian will be derived from Fermi's golden rule and has the 
following form 9 



p r +co 

1{ y) = J dte^ h ([H T (t),H, 



+ 00 
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(18) 



The fact that the magnetophonon mode is now gapped, 
contrary to the single layer case in which is gapless, is an 
important distinction between the single layer and bilayer 
case for the tunneling characteristics. 



where the correlation function definitions are 

I=F(i) = (H T (t)H+(0)), 
I^it) = (H+(t)H T (0)), 



(19) 
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and the time-dependence is meant in the interaction pic- 
ture representation, namely 



(20) 



For the calculation of the correlation functions in Eq. 
()19|) we use the same approach as JK. We assume the 
tunneling process is statistically independent from the 
collective mode propagation and can thus be averaged 
independently. For the statistical averaging of the collec- 
tive modes we use the linked cluster expansion method^ 
which turns out to involve only one link in the expo- 
nential resummation. In Appendix D we show in more 
detail how the calculation proceeds. For the correlation 
function we obtain 



lT(t) = v(\-v)T 2 C{t) 



(21) 



where 
C(t)=exp{-£ 
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JQ.t 
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and N s is the boson thermal occupation number for the 
collective modes of the system. To get the form of Eq. 
(jT§|) it suffices to interchange A and B in Eq. No- 
tice that the in-phase modes drop out from the above 
correlation function due to their matrix element prop- 
erty M s a — M s b, this is not the case for the out-of- 
phase modes which, as we show in the Appendix, obey 

M sA = -MsB- 

The experimental temperature range in the tunneling 
current results is of the order of 0.1 K^10~ 5 eV while the 
bias voltage is in the range of mV, so a zero temperature 
calculation is appropriate; this simplifies things consid- 
erably since the bosonic occupation numbers N s = 0. In 
addition, due to the high magnetic field the magnetoplas- 
mon modes will have a large gap and will not contribute 
to the electron coupling. Gathering all this together and 
switching to dimensionless units for the momentum in- 
tegration (x = q/qo) we find the following result for the 
time dependent correlation function 



C(t) = cxp 



dxf(x)( 



3 —ioj(x)t 



i 



(23) 



In Appendix D we provide the analytic expression for the 
weight-function f(x), while ui(x) is given by Eq. (|13fl or 
specifically by Eq. (|D6|) . By differentiating Eq. 
and taking the Fourier transform we obtain the follow- 
ing equation for the Fourier transform of the correlation 
function 



u)C(w) = I dxf(x)oj(x)C(ui — ui(x)). 



(24) 



As we show in Appendix D this correlation function is 
zero for uj < 0. 



The integral equation in Eq. (|24(l is difficult to solve 
analytically. However, we can find the asymptotic be- 
havior of C(ui) by expanding in lu(x). To lowest order we 
obtain a first order differential equation with the solution 



C(w) 



exp 



(u - ci) 



2 1 



2c 2 



(25) 



where we provide analytic expressions for c\ and ci in 
Appendix D. This expansion gives us an analytic pre- 
diction for the peak bias voltage, namely c±, and the 
general behavior of the correlation function at large bias 
voltages. We can investigate the behavior of c\ in the 
two limiting cases when the interlayer Coulomb inter- 
action is much weaker than the intralayer one and vice 
versa. We find that for the d ^> a,Q limit, c% scales as 
1/ao, and the width of the Coulomb barrier peak ^Jcq, 



scales as ^foi ~ 1/a^vB. This is the same phenomeno- 
logical behavior JK extract from their calculations. In 
the opposite limit d <C do, we find that c\ ~ d 2 ja\ and 
y/7-2 ~ d/flQ-v/S. This limit is absent in the JK model. 
We should notice here that in this regime the actual be- 
havior of the system will be significantly modified by co- 
herence effects but those effects will be pronounced at 
zero-bias. At finite bias values we see from experimental 
results^ that the Coulomb barrier peak survives, but it 
is "red-shifted" significantly which, as we see, our model 
is able to reproduce as a limiting behavior. 

The importance of the asymptotic result lies on the 
insight it provides on the tunneling current behavior. We 
use this information to build an Ansatz for the correlation 
function. We are after a qualitative analytic expression 
that will capture the basic physics and reproduce the 
Gaussian asymptotic behavior. We choose the following 
Ansatz 



C(w) = Nu r e 



-Au 2 



(26) 



As we show in Appendix E, this choice reproduces very 
accurately the numerical solution of the equation which 
is our strongest justification for using it. For the pa- 
rameters involved in this choice of C(ui) we develop a 
self consistent method of evaluation. This method is 
based on the moment expansion of C{ui). If we multiply- 
differentiate Eq. ll'-'il) we end up to the following moment 
equations 



/>oo 

/ dujC(uj) = 2tt, 
Jo 

du>LoC(uj) = 2nci, 



duju 2 C(uj) = 2tt(c 2 + cf). 



(27) 
(28) 
(29) 



These are the three equations that our Ansatz should 
obey self-consistently. We show in Appendix E in more 
detail how we extract the values for N, r and A from 
them. Our theoretical value for the peak bias voltage 
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can be easily calculated now from Eq. I|26|) . The result 
is 



V = 



lOOOfi 



(30) 



where we have converted it to mV. In Fig. Q we present 
the results the moment expansion method produces. 
The theoretical parameters are directly taken from the 
experiment^ The bilayer sample area is S=0. 0625mm 2 
and the single layer electron density is no =1.6 x 10 11 
cm~ 2 . The perpendicular magnetic field varies from 8 
T to 13.75 T and the Wigner crystal lattice parame- 
ter has the value ciq ~ 270 A which corresponds to the 
stable hexagonal lattice configuration (n^ = 2/\fia^) 15 
The double well separation distance is <i=175 A and the 
dielectric constant of GaAs is e=12.9eo ~1.14 x 10~ n 
F/m. The transverse sound velocity for the electron gas 
is given byi^ 



c T 



'0.0363 



V3. 



emao 



53552 m/s. 



(31) 



For the electron mass we use the electron effective mass 
value in the GaAs background m=0.067m e . For the K 
parameter we notice that by construction K/uq is a mea- 
sure of the energy density per electron due to the short 
range part of the interlayer Coulomb interactions. So we 
assume it will scale accordingly as 



K 
n 



2 1 Aired 



K > 0. 



(32) 



The value of k is a measure of the magnetophonon gap 
in this system. Cote and collaborators have performed 
time-dependent Hartree Fock calculations investigating 
the magnetophonon dispersion relation and crystalline 
phases for the bilayer system^ They were able to pro- 
vide us with a k=0.0085 value. With this last piece 
of the puzzle in place we are able to produce our ana- 
lytic results for the tunneling current shown in Fig. ffl. 
As we see, the model captures the qualitative behavior 
of the current and it accurately reproduces the experi- 
mental results for the peak bias voltages associated with 
different applied magnetic field values. Notice though, 
that the pronounced experimental low bias suppression 7 
is lifted in our case although present in the JK model 9 
and in spectral weight calculations 17 associated with this 
bilayer system. We attribute this behavior to the very 
small value of k and mostly to the role of interlayer in- 
teractions. In particular, as it can be scon in Fig. (JTJ, the 
Coulomb blockade tends to be restored with increasing 
magnetic field because the role of interlayer interactions 
is diminished in that limit. We observe the same behav- 
ior as well when the interlayer distance d is increased, 
which signifies the same physical limit, where essentially 
the JK results are recovered. We believe that in order to 
capture all quantitative characteristics of the tunneling 
current for the whole range of bias voltages one has to go 
beyond the continuum approximation employed here. 
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FIG. 1: Tunneling current curves for different magnetic field 
values using the moment expansion solution of Eq. fZQ. The 
legend shows the peak bias values calculated by Eq. I3UII . 
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FIG. 2: Normalized tunneling current curves for different in- 
terlayer separation distances d measured in A. 



Another way of calculating the correlation function is 
by numerically integrating Eq. I|24(l . In Appendix E 
we report analytically on this approach. The results we 
get are shown in Fig. J3J and are very similar to the 
moment expansion results; the same applies for the peak 
bias values. 

At this point we are in position to analytically inves- 
tigate how our solution behaves for different interlayer 
separation distances d. Notice that the dependence in d 
is introduced in two places in our model, the one is the 
tunneling matrix elements T (which depends exponen- 
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tially on d) and the other is through the long-range part 
of the interlayer Coulomb interaction. Since we are inter- 
ested in the latter only, we will normalize the I-V curves 
for different d values to the peak current value (this has 
the effect of removing the tunneling matrix element con- 
tribution). According to experimental results^ we expect 
to see a "red-shift" in the I-V curves associated with the 
bounding interlayer exciton energy —e 2 fed. Our results 
are shown in Fig. J5J and as we see our model captures 
this kind of important interlayer correlation physics. Ac- 
cording to our model, we can conclude that the effect 
of the interlayer Coulomb interaction in the real system 
is two-fold. First, the short-range correlations derived 
from the interlayer interaction introduce a gap in the low 
frequency excitations, which is well pronounced in the 
experimental results,- and second, the long-range effects 
of the interlayer interactions "soften" the response of the 
bilayer system to tunneling, through the exciton creation 
associated with tunneling events. 



V. SUMMARY 

We have introduced a model for tunneling in a bilayer 
system in high magnetic fields, which explicitly includes 
the interlayer interactions. The bilayer system of elec- 
trons is modeled as two correlated Wigner crystals in the 
continuum limit. The low energy collective modes of the 
crystals are coupled to the tunneling electrons in the sys- 
tem, providing an energy dissipation mechanism to relax 
the charge defect produced by the tunneling event. This 
model allows for an approximate analytic solution to the 
tunneling current problem that captures the basic experi- 
mental behavior of the tunneling current curve and fairly 
accurately reproduces the peak bias values for different 
applied magnetic field cases along with the expected I-V 
dependence with interlayer distance d. In the future we 
would like to find a way to include coherence into the 
model and reproduce the zero bias peak anomaly^ the 
tunneling current exhibits. 
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APPENDIX A: SINGLE LAYER EIGENMODES 

We have expressed the single layer dynamics by writing 
down the Lagrangian of Eq. J2J). The canonical quantiza- 
tion process is a textbook exercise where one has to solve 



the Lagrange equation of motion for the fields, namely 



ut + oj c ul + ut u t = 0, 
u L - u c u T + lo 2 l u l = 0. 



(Al) 
(A2) 



Notice the effect that the magnetic field has in coupling 
the otherwise orthogonal longitudinal and transverse zero 
field eigenmodes. The above system of equations is easy 
to prove that has eigenvalues given by Eq. JSJ. The 
general solution for the eigenmodes is of the form 



Ut 



= A T e iu + t +A're- iu '+ t +B T e iu '- t +B' 1 



[-A T ^ w+i + A' T e- iu '+ t } 



+ i- 



-B T e u 



(A3) 



(A4) 



where At, A' t and Bt, B' t are the coefficients asso- 
ciated with the creation and annihilation of the two 
modes respectively. These are directly related to the 
creation and annihilation operators by the normaliza- 
tion condition. The way we reach to the latter is by 
expressing all the higher derivative terms in terms of the 
canonical displacements and the canonical momenta at 
time equal to zero which obey the commutation rela- 
tions [ui (q),pj (<i')] — ih8ij{2ir) 2 5 2 ((\ — q'). After a 
fair amount of algebra we have convinced ourselves that 
the operators obey the correct eigenmode commutation 
relation. By applying the normalization condition we can 
reach the final form of the eigenmode operators given by 
(implicit q dependence will be suppressed unless neces- 
sary) 



ai(q) 



«i 



+ co+l cu 



nQin ' 
1 



-PL 



-UL - 



+ 



TIqTYL 



T> UT ° 



-PT 



(A5) 
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where they obey [aj(q), aj(q')] = (2-7r) 2 <5 2 (q— q')5y and 
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We can finally express the displacement fields in terms 
of the eigenmodes in Schrodinger's representation as 



u T (q) = 




k {a\+ ai ) 



(A9) 
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APPENDIX B: BILAYER EIGENMODES 

Starting with Eq. JSJ) and introducing the in-phase v = 
(ug + ua)/2 and out-of-phase u = ub — u.a modes, we 
find that the Lagrangian can be written as L — L m +L ou t, 
with 



Lm = no 



d 2 q 1 



(2tt) 2 2 

2 2 -o2 2 



-(2m) v\ + i'T + lo c (v t v l - v l vt) 



Lout — 



where 



d 2 q 1 m 



(2tt) 2 2 2 



2 

rp (J,rp 



Ql = c 2 L q 2 



2K 



mriQ 



ol 



2 2 

Cl<7 



(Bl) 

+ + IjJ c (utUl — ULUt) 

(B2) 

(B3) 

(B4) 

(B5) 
(B6) 



2 2 



2Jf 



rrmo 



e 2 n 



2me 



9(1 -e 



e 2 n 



2me 

Oj 1 = 



9(1 



We see that the four modes have decoupled into two in- 
phase (with total mass 2m) and two out-of-phase (with 
reduced mass m/2) modes that are governed by effec- 
tive single layer Lagrangian dynamics. Notice that the 
in-phase transverse mode remained unchanged. In order 
to obtain the analytic results for the out-of-phase modes 
displacement field operators we have to redefine the pa- 
rameters in Eqs. (|A5IA10|) as m — > m/2, (jj, — > f2y, 
lu l — > £l 2 L while for the in-phase mode we have m — » 2m, 
w| — > 0\ . The new eigenvalues f2± will be given by Eq. 

and a corresponding one from Eq. (J5J for 0± when we 
carry out the above frequency changes. Notice that the 
mass in the cyclotron frequency formula does not change 
in either case. We will use the operator form of these 
eigenmodes to calculate below their coupling matrix ele- 
ments with the electrons in the system. 



APPENDIX C: 



COLLECTIVE MODE COUPLING 
MATRIX 



The coupling between a point charge with charge den- 
sity n e (r) and a charge fluctuation Sn(r) in either of the 



two layers is given by 



** **** r _ r . 



n e (r)6n A (r A ) 



e 

Ant 



+ — ld 2 r ld 2 r B 



n e (r)5n B (vB) 



(CI) 



In the continuum approximation to lowest order we have 
Sua = — JioV ■ Uyi, and we can place the electron at the 
origin so that n e (r) = <5^ 2 ^(r). After introducing the in- 
phase and out-of-phase displacement fields and Fourier 
transforming, the coupling becomes 

H cpl = (c A c A -J B c B ) ^(l-e-^iq.u, (C2) 

where we have dropped a constant term (which depends 
upon v). Next, we use the operator definition for the 
out-of-phase displacement field found previously 



u L = -i/ioj + ifiai - if 2 a\ + «/2«2, 



where 



/i = 




n m (Q 2 + -^ll)(n 2 + -Q 2 _) , 



0_ 



n m (n 2 L - n 2 _){n 2 + - n 2 



(C3) 

(C4) 
(C5) 



and the coupling term assumes the form of Eq. I|15|l . 
where the matrix elements are given by 



MtA = -M. B = -^(l-e-"" 



)fs, (C6) 



and s = 1, 2 refers to the two out-of-phase bilayer eigen- 
modes but any summation on it is meant to include an 
integration in q as well. 



APPENDIX D: CORRELATION FUNCTION 

To calculate the correlation function given by Eq. 119|) 
we use the linked cluster expansion method^ As it turns 
out there is only one link involved in the exponential re- 
summation of the series. As pointed out by JK, if we 
assume that the tunneling events are faster than the re- 
laxation time for WCs, then the thermal averaging of the 
magnetomodes and the tunneling electrons is statistically 
independent. The correlation function becomes 



F(t)=T 2 v{l-v){T t ex V 



(H f - Hi 



- / dt'e— 



(Dl) 
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where 

Hi = e A -A A - 



H 



f 



E 



f£l. 



iM sA 
Ml s 

iM sB 
Ml, 



iMsB 



), (D2) 



(D3) 



are the Hamiltonians of the system before and after a 
tunneling event respectively. Using the linked cluster ex- 
pansion method we obtain the result of Eq. (|21|) . aside 
from an overall phase factor (e^ — A^— e^+A^i/ft which 
we assume to be zero. As JK point out, this amounts to a 
complete Wigner crystal relaxation after each tunneling 
event, where the polaron shifts A^ — A^ defined by 



E 



\M aA{B) \ 



(D4) 



The parameter a is dimensionless and gives a measure 
of the magnetophonon gap. We have used Eq. lE!^ to 
produce the final equation that defines it. The parameter 
7 is dimensionless and measures the relative interaction 
strength between the interlayer and intralayer Coulomb 
interactions. The parameter (3 is dimensionless as well 
and does not depend on the Wigner crystal lattice pa- 
rameter ao as it turns out if the dependence of ct given 
by Eq. (|31|l is taken into account. In order to derive 
the result for the asymptotic behavior of the correlation 
function we have to Taylor-expand C(lu) involved in the 
integral of Eq. (|24[1 in powers of ui(x). It is straightfor- 
ward to show that to lowest order the first order differen- 
tial equation we obtain has the Gaussian solution given 
by Eq. 1|25|) . The coefficients involved in the expansion 
are given by 



compensate the Madelung energies e B — of the inter- 
stitial positions. Another point of interest is that as we 
see from Eq. (|22[) the matrix elements of the electron 
coupling with the collective modes of the system enter in 
the correlation function formula as a difference, so that 
the tunneling electron only couples to the out-of-phase 
modes. When an electron tunnels from one layer to the 
other the remaining electrons in that layer have to fill in 
the hole; at the same time the electrons in the receiv- 
ing layer need to "open up" and create a hole for the 
incoming electron, so we see that the relative motion of 
the two charge densities should be out-of-phase during a 
tunneling event. 

Since the magnetoplasmons have a large gap in high 
fields we only need to consider the magnetophonon ma- 
trix clement, which in the high magnetic field case will 
be given by 



/(*) 



x(l 



-7Z\2 



1 



4<?§ S + x 2 + 2a + 0x(l - e-f") ^ 



1 



1 



[a + /3x(l - e-~< x )} 3 / 2 S-a-x 2 ' 
For the magnetophonon frequency we will have 



L0{x) 



2 2 



yja + f3x(l -e—t x )\fx 2 



(D5) 



(D6) 



For convenience we have introduced dimensionless units 
for the momenta (x — q/qo), and we have set cl = 0. 
The parameters in the equations above are 



n 
8nhm 



2K 



mnoc^qg 2ir 2 mc 2 ^ {Iqo) 2 ed' 



= 



e 2 n 



5 = 



7 = dq , 

( UJc 
\ct9o 



(D7) 

(D8) 

(D9) 
(D10) 
(Dll) 



ci = / dxf(x)u>(x) 
Jo 



,4 pi 



dx- 



c(l - e-T x ) 2 



?t% Jo 5 + x 2 + 2a + (3x(l - e-~< x ) 



a + f3x(l — e ~ ,x ) 8 — a — x 2 



(D12) 



C2 = / dxf(x)u) (x) 
Jo 



3 r 1 



= c- 



dx- 



c(l - e"T x ) 2 



4#> Jo 5 + x 2 + 2a + (3x(l - e-~< x ) 



x z + a 



1 



a + (3x{l — e ryx ) 5 — a — x 



2 ' 



(D13) 



As we mentioned in the text we can expand to one or- 
der higher and analytically solve the differential equation 
as well, but the correction will not be very useful since 
these results are valid only in the asymptotic region that 
a lowest order calculation suffices to catch the essential 
characteristics. In the d 3> ao case we can ignore the ex- 
ponentials in the integrand of Eqs. i|D12l |D13|) and as it 
turns out a c\ ~ 1/ao and ~ l/a^B scaling behavior 
emerges. In the opposite limit d <ti ao, if we expand the 
exponentials in the integrand of Eqs. (|D12llD13|l we find 
that ci ~ d 2 / a o an d \fc2 ~ d/a^B. 

At this point we would like to justify the assumption 
that the correlation function is zero for zero or negative 
bias voltage (uj < 0). By expanding the exponential in 
Eq. (|24|l and interchanging summation with integration 
we can perform the last integral that gives a delta func- 
tion in the frequencies. Since the magnetophonons are 
always positive definite (gapped) we see that the corre- 
lation function 

oo „i „i 
C(w) =27re-/o *«/(«) E ~ / ^i/^i) ■ • • / dx n f(x n ) 

n=0 H ' ^° J° 



x5(oj — w(xi) 



■ CO 



(x n )) 



(D14) 
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has to be zero for any uj less than the lowest possible 
magnetophonon frequency value. 

With this last piece of information about the corre- 
lation function we can proceed now into developing our 
method of moment expansion. As we have already justi- 
fied we use the Ansatz of Eq. H26JI and the three different 
moments of the correlation function given by Eqs. (|27h 
I29|l . At this point it is convenient to make a change of 
variable ui = eVyiOOOft to convert the argument of C(oj) 
into mV and define the corresponding parameter in the 

exponent as A = A( 10 Q 0rt ) ■ By performing the integrals 
the three moment equations become 



1 



-N 



1000ft 



2 



1000ft 



1000ft 

r+3 

A 



r+l 

A" 

r+2 



1 



2tt, 



= 2nci, 



= 2tt(c 2 



(D15) 



(D16) 



(D17) 



We can divide Eq. l|D16)l and Eq. (|PT7)l by Eq. l|D15l) 

and then equating the two we reach the self consistent 
equation for r, namely 



r(f) 



1 + c 2 /c{ 



(D18) 



We numerically solve the above equation and obtain a 
value for r. The usual range of r for the magnetic field 
values considered is 1/2 < r < 2. This can be regarded as 
an estimate for the low-bias current power-low behavior. 
Having r at hand we can go back and evaluate the rest of 
the parameters. The final form of the correlation function 
becomes 



C(V) = N 



1000ft 



V r e 



-AV* 



(D19) 



where V measures in mV. To obtain the peak bias value 
result of Eq. (|30[) we find the root of the first derivative 
of the above equation. 



APPENDIX E: NUMERICAL INTEGRATION 

We have numerically integrated the integral equation 
for the correlation function in two different ways — first 
by a direct integration of the integral equation and then 
by introducing the density of states (similar to the JK 
method) . As both methods give the same results we will 
only report on the latter. The equation to be integrated 
is 



zC(z) 



dxf(x)^C(z- ^ 



Jo 



7o 



uj(x) -, 
7o ' 



(El) 



lOOOfi 



to 



where we have introduced the parameter 70 
convert the frequency argument of the correlation func- 
tion into mV. Before we proceed we should notice that 



the magnetophonon frequency is bounded in a region 
c^i < < a 2 which means that the density of states 
is non zero only in that range. The values of ot\ and a 2 
are given by substituting x = and x — 1 into Eq. I|D6|) 
respectively. The upper bound a 2 appears due to the 
momentum cutoff we have introduced. The final form 
for the correlation function integral equation will be 



( 1 r z 



C(z) = { 



J ai d y9(y)C( z - V), ai<z<a 2 , 



(E2) 



^J a * d y9{y)C{z-y), z>a 2 , 



where the definition for the density of states is 



Ay) = y 



1 I dm(x) 
70 I dx 



(E3) 



and x{y) is the root of the equation u>(x) = "foy. In this 
approach one has to "jump-start" the algorithm with an 
assumption for the low bias points. We can show that 
for small frequencies (but greater than the lower mag- 
netophonon bound) the correlation function will have a 
z — a.\ behavior. In Fig. © we show our results for the 
tunneling current using this method. As we see the quali- 
tative behavior of the tunneling current is very similar to 
our analytic result, and the peak bias values are similar 
as well. 



I 

'*. 13.75 T 




-Peak = 


5.59 mV 


-■Peak = 


6.32 mV 


■- Peak = 


6.84 mV 


■■■ Peak = 


7.52 mV 



10 15 20 25 

Interlayer Voltage (mV) 

FIG. 3: Tunneling current curves for different magnetic field 
values produced by numerically integrating Eq. 1241 . The leg- 
end shows the peak bias values obtained with this approach. 
They are in close agreement with the analytic results. 
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